How branching can change the conductance of balhstic semiconductor devices 
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We demonstrate that branching of the electron flow in semiconductor nanostructures can strongly 
affect macroscopic transport quantities and can significantly change their dependence on external 
parameters compared to the ideal ballistic case even when the system size is much smaller than 
the mean free path. In a corner-shaped ballistic device based on a GaAs/AlGaAs two-dimensional 
electron gas we observe a splitting of the commensurability peaks in the magnetoresistance curve. 
We show that a model which includes a random disorder potential of the two-dimensional electron 
gas can account for the random splitting of the peaks that result from the collimation of the electron 
beam. The shape of the splitting depends on the particular realization of the disorder potential. At 
the same time magnetic focusing peaks are largely unaffected by the disorder potential. 
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I. INTRODUCTION 



Transport of two-dimensional electrons in state-of-the- 
art modulation-doped semiconductor heterostructures 
suffers mainly from small angle scattering off charge fluc- 
tuations in the donor layer, which is separated from the 
two-dimensional electron system by a spacer layer——. 
The amplitude of the disorder is only a few percent 
of the Fermi energy, and the electron mean free path 
can exceed hundreds of microns^^^. Under these circum- 
stances, transport through nanostructures much smaller 
than the mean free path is assumed ballistic, i.e. in- 
dividual electrons follow almost perfectly the paths pre- 
scribed by Newton's law under the external applied forces 
and the forces associated with the nanostructure's con- 
finement potential. Impurity scattering is thought of as 
negligible. Nevertheless, experiments, which used spa- 
tially resolved recordings of the change in conductance 
induced by a charged tip of a scanning probe microscope, 
showed that the current density emerging from a quan- 
tum point contact was branched as a result of small angle 
scattering onl}^ — . The reported conductance changes in 
these and further experiments using the same measure- 
ment techniques were much smaller than the conductance 
quantum e^/h. At the same time many experiments 
have been interpreted in terms of purely ballistic effects, 
e.g. focusing by electrostatic^ and magnetic^^ fields or 
wall geometric^ even in semiconductor materials with 
shorter mean free paths. It could thus seem that branch- 
ing has little impact on macroscopic transport quantities 
and hence is of no relevance for most transport experi- 
ments. 

In this article, however, we use a novel magnetic fo- 
cusing device to demonstrate both experimentally and 
theoretically that the macroscopic transport quantities 




FIG. 1. The first magnetoresistance peaks recorded for two 
mesoscopic devices fabricated from the same heterostructure. 
The inset shows the scanning electron microscope picture of 
the novel magnetic focusing device. Bq is the magnetic field 
at which the cyclotron radius is commensurate with the wall 
length of the device. 



of nanostructures can in fact be strongly influenced by 
branching. This poses the question of how ballistic trans- 
port effects observed in many supposedly ballistic exper- 
iments actually are. To illustrate this let us examine the 
two sections of magnetoresistance curves from two differ- 
ent samples fabricated to the exact same specifications 
shown in Fig. [T] (details will be given below). In an ideal 
(i.e. ballistic) sample patterned in this specific layout 
(depicted in the inset) electrons emitted from point con- 
tact Ci are deflected by a magnetic field and for a certain 
value B = Bq will be directed towards the second point 
contact C2. Therefore, one expects to observe a single 
peak in an appropriate transport quantity at B/Bq = 1. 
Sample 1 shows exactly this behavior, sample 2, however, 
shows an unexpected splitting of the peak. So far, many 
experimentalists would argue that sample 2 is a defective 
sample, where an unfavorably located impurity spoils the 
measurement and thus that the sample should be dis- 




FIG. 2. Branching of a plane wave front propagating in x- 
direction. (a) The electron wave flow intensity approximated 
by classical ray dynamics is plotted in gray scale as the plane 
wave propagates along the x-direction in the disorder land- 
scape shown as a green/white color rendition in the back- 
ground. The force in x-direction is ignored (quasi-2D). The 
wave front remains a vertical line in coordinate space. Verti- 
cal lines mark the position of the plane wave front at times 
ti, t2, . . .,t9. The bottom panel displays the wave fronts at 
these times in phase space {y,Vy). Caustics are identified at 
the turning points (purple dots), (b) Flow density p as a 
function of y at times and tr (vertical cross-sections of the 
flow density plot in a). Also shown at the bottom are the 
corresponding wave fronts in phase space. The flow density 
peaks at the caustics. 

carded. We show, however, that curve 1 and 2 are both 
fully compatible with exactly the same amount of impu- 
rity scattering and that it can not be argued that sample 
1 is in any way better than sample 2. On the contrary, 
even though both samples are actually extremely clean 
devices, it is a rather fortunate coincidence that curve 
1 agrees well with the expectations for an ideal system. 
In addition, our new magnetic focusing device allows us 
to explain why these consequences of branching have not 
been seen in previous magnetic focusing geometries. 



II. CAUSTIC FORMATION BY DISORDER 
AND A MAGNETIC FIELD 

A. Branching of an initially plane wave 

We prelude the detailed presentation of our results by 
an illustrative explanation of the mechanism of branch 
formation, as depicted in Fig. [2^ and b. An initially ho- 
mogeneous particle flow or plane wave front (restricted 
source in momentum space) is propagated along the x- 
direction through a two-dimensional disorder potential. 
A color rendition (from green to white) of the disorder 
potential is plotted in the background of panel a. For 
the sake of simplicity, only the electrostatic force exerted 
by the disorder potential in the ^/-direction is considered. 
The force in x-direction is ignored, so that the longitudi- 
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FIG. 3. Influence of disorder on the flow density p (gray scale 
rendition) emitted from a point source, (a) gray scale rendi- 
tion of the flow density p in the (x, ^)-plane when particles are 
emitted from a point source with a cosinusoidal angular dis- 
tribution. They are subjected to the disorder potential color- 
coded in green and white. Forces in both x and y-direction 
were taken into account and the magnetic field is absent, (b) 
The formation of caustics in a transverse magnetic focusing 
geometry. Particles are emitted from the top point contact. 
The flow density is plotted on a gray scale at a magnetic field 
where a caustic forms at the collecting point contact, (c) 
Collimation (left panel) at B = Bq and magnetic focusing at 
B = 2Bo (right panel) in a corner device consisting of two 
quantum point contacts placed at a 90° angle. Red curves 
in the lower panels show the flow density hitting the lower 
wall, (d) Same calculations as in c but in the presence of 
weak disorder. The standard deviation of the amplitude of 
the disorder potential corresponds to 2%Ef. 



nal velocity Vx stays constant and the wave front remains 
a vertical line in coordinate space. This model makes the 
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principles of branching particularly easy to understand 
and yet captures all of its important featured. Later 
on, an extension of the model will also allow us to study 
the statistics of the formation of branches in a magnetic 
field analytically. The particle or fiow density p is shown 
in panel [2^ using a gray scale. While p is initially ho- 
mogeneous or independent of it develops features at 
later times as illustrated in panel b where p{ii) is plot- 
ted at selected times. At a time between ^4 and ^5, a 
strong peak develops in p. This heralds the first branch. 
In the phase space {y^Vy)^ the wave front develops a pair 
of initially coalescing turning points, which subsequently 
separate (panel a bottom). Here, the classical ray density 
in coordinate space diverges and also the quantum me- 
chanical wave intensity would peak nearby. In between 
the turning points, the wave front folds in phase space 
and covers the same spatial coordinate three times. As 
a result, the local density is enhanced. The path traced 
by a turning point constitutes a fold line or caustic. A 
branch is referred to as the spatial region in between two 
such random caustics. 



B. Branching and magnetic focusing with a point 
source 

Figure [3^1 illustrates another instance of caustic for- 
mation and branching when particles are emitted from 
a point source (restricted source in coordinate space). 
In this example, no approximation concerning the elec- 
trostatic force associated with the disorder potential is 
made and the full two-dimensional particle dynamics is 
considered. Branches appear on similar length scales as 
for the simplified plane wave case due to the same basic 
mechanism^^. This setup can be implemented straight- 
forwardly in a GaAs based 2DEG on which a quantum 
point contact (QPC) is patterned either by etching or the 
split gate technique^ —. 

It is instructive to oppose the appearance of branches 
due to the disorder induced formation of random caus- 
tics to the focusing of two-dimensional electrons emitted 
from a point source in a perpendicular magnetic field B 
in the absence of disorder. We do this, because also in 
the magnetic focusing problem caustics play an impor- 
tant role^iiSi^S. In a magnetic field, the electrons execute 
circular cyclotron orbits with a radius r = ^^^^^ , where 
n is the electron density. When emitted from a point 
source, the electron trajectories converge at a distance of 
one cyclotron diameter away from the point source and 
a caustic forms. The enhanced local current density can 
be detected with the help of a second collecting quantum 
point contact in the transverse magnetic focusing geom- 
etry shown in Fig. ^pi^^^. The collecting QPC is at a 
distance a from the emitting contact. The enhanced lo- 
cal density due to magnetic focusing arises at a magnetic 
field for which the cyclotron diameter 2r equals a. When 
enforcing zero net current flow through this QPC and 
measuring the voltage drop across, the enhanced local 



density can be detected as a voltage or resistance peakiS. 
One may anticipate that disorder induced branching also 
produces resistance peaks when sweeping the magnetic 
field. 



C. Focusing in the corner device 

In order to search for evidence of branching in dc bal- 
listic transport we have chosen the corner shaped device 
depicted in Fig. [3fc, which, as we will show, can distin- 
guish between deterministic magnetic focusing at the col- 
lecting QPC and random focusing caused by branching. 
The left panel shows trajectories in the absence of dis- 
order at the field Bq for which r = a. The fiow density 
along the bottom boundary is plotted and forms a broad 
peak centered around the collecting point contact open- 
ing for this magnetic field. It is not caused by magnetic 
focusing since a caustic has not yet developed. Rather 
it results from the collimating properties of the emitting 
point contact from which electron trajectories leave with 
a cosinusoidal angular distribution. One may anticipate 
that it produces a first resistance peakiii2^. At the field 
2Bo (right panel of Figure [Sb), the particle fiow converges 
first into the device corner where a caustic has formed 
due to magnetic focusing. The corner was chamfered to 
have well defined refiection. After specular refiection the 
electron trajectories refocus at the collecting QPC giving 
rise to another resistance maximum. 

Figures [3b and d illustrate how the disorder potential 
affects the particle fiow for both values of the magnetic 
field. While at 2Bo the infiuence of disorder is weak, the 
fiow density has been drastically altered for B = Bq and 
now exhibits multiple maxima of comparable size. The 
magnetic focusing feature at 2Bo remains largely unal- 
tered and is more robust against disorder induced branch- 
ing. These conclusions can be generalized and also hold 
at higher magnetic fields. In the absence of disorder, the 
fiow density reaches a maximum at the collecting QPC 
for B = kBo, where A:=l,2,. . .. Collimation and focusing 
are responsible for these maxima at odd and even values 
of k respectively. The collimation features are prone to 
disorder induced splitting, while the focusing features are 
generally more resilient. As opposed to the conventional 
transverse magnetic focusing geometryiS, in this corner 
device deterministic focusing and collimation are sepa- 
rated on the magnetic field axis and hence this geometry 
lends itself particularly well to confirm the theoretical 
predictions. 



III. EXPERIMENTAL RESULTS 

The devices are fabricated from a modulation doped 
GaAs/AlGaAs heterostructure in which the 2DEG is lo- 
cated 150 nm underneath the crystal surface. Split gates, 
arranged as shown in the insets of Fig. HI form a 90° cor- 
ner with a QPC along each leg. One split gate is shared 
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FIG. 4. Magnetoresistance traces measured in corner devices with different QPC resistances, Rqpc- The electron mean free 
path of the 2DEG is 45 /xm. (a-b) Magnetoresistance traces for two devices with a chamfered corner fabricated from the same 
wafer show different behavior of the resistance i^focusing feature near B/Bq—1. In each panel, curves are shown for various 
values of the resistance of the quantum point contact, i?QPc-(c) Magnetoresistance data for a device with a sharp corner design. 



by both QPCs and defines the corner. Devices with a 
chamfered corner (Fig. [4^ and b) as well as a straight 
corner (Fig. |4] c) are investigated. The electron density 
n equals 2.5-10^^cm~^ and 2.2-10^^cm~^ in the devices 
used in Fig. |4^,b, and c respectively. The wall length a 
from each QPC to the corner is 3 /im. The electron mean 
free path of the 2DEG is 45 /im, one order of magnitude 
larger than the ballistic electron trajectories relevant for 
these studies. Transport measurements in a perpendicu- 
lar magnetic field are carried out at 1.4 K by driving a 
sinusoidal current / of 5 nA with a frequency of 13.3 Hz 
through the injector QPC. The gate voltage applied to 
all three gates defining the QPCs is identical. The volt- 
age drop across the detector QPC, V ^ and the injector 
QPC are detected with separate voltage probes using a 
lock-in technique (measurement configuration shown in 
the insets of Fig. |4]). 

Three typical experimental data sets are plotted in 
Fig. |4j The resistance data i^focusing = V/I are recorded 
on devices which have an identical size but possess dif- 
ferent realizations of the disorder. Since all devices are 
fabricated from the same heterostructure, the disorder is 
characterized by approximately the same statistical pa- 
rameters. The curves within each panel are acquired for 
different gate voltage, i.e. different values of the QPC 
resistances. For the bottom curves at low QPC resis- 
tances many modes propagat o^^i^^ . The top curves are 
recorded for transmission of a single mode. The colli- 
mation features at odd values of B/Bq are in general 
much broader than the focusing features at even values 
in agreement with the broad distribution for the particle 
flow density in the left panel of Fig. [3t. They frequently 



split into two peaks depending on the device, i.e. the 
specific disorder realization (for instance at B / Bq=1 in 
panel b and at B / Bo=3 in panel c). We assert that the 
splitting of the collimation features is a result of disorder 
induced branching as in the left panel of Fig. [3t. Such 
splitting is absent for the focusing features at even values 
of B/Bq. The focusing features are more robust against 
variations of the disorder potential. The data in panel c 
were recorded on a device with a corner that is not cham- 
fered. Processing as well as depletion will smoothen this 
corner somewhat however the direction of the specular re- 
flection is not as well defined and focusing features which 
require specular reflection in this corner have dropped in 
amplitude. 



IV. TRANSPORT SIMULATION IN THE 
CORNER DEVICE 

To corroborate our assertion that the splitting of the 
collimation features comes from disorder induced branch- 
ing, we have numerically calculated resistances for differ- 
ent disorder potential realizations. We simulate trans- 
port in the corner device by following classical trajecto- 
ries from the emitting QPC until they either reach the 
collecting QPC and contribute to the transmission, or 
until they leave the system to the right of the collec- 
tor. Within the Landauer-Biittiker- formalism, we ap- 
proximate the resistance i^focusing measured in a four- 
terminal setup to be proportional to the transmission 
probability from emitter to collector. The QPC is tai- 
lored as a sum of variable-depth hyperbolic tangents, and 
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FIG. 5. Numerical transport simulation and its compari- 
son with experimental results, (a) The number of trajectories 
which reach the collecting QPC as a function of B / Bq in dif- 
ferent realizations of the disorder potential in a corner device 
with chamfered corner. Peaks at odd multiples of Bq are 
strongly affected by branching. The inset shows a contour 
plot of the potential used to simulate the QPC. (b-c) Com- 
parison of experimental data (magnetotransport curves for 
Rqpc ~ 3 kQ in Figs.[3]3 andfH:.) and numerical simulations 
for different disorder realizations. 



we assume the presence of a saddle potential inside the 
QPC. Soft wall effects due to depletion are modeled by 
using a quadratic potent ial^S. In the simulations, it is as- 
sumed that the particles enter the QPC from a lead with 
a cosinusoidal angular distribution. The saddle potential 
then has the effect of collimating the flow. We find that 
to obtain a collimation peak which is consistent with the 
experimental results, a saddle potential of approximately 
SO%Ef is needed. The precise functional form of the sad- 
dle potential is given in [26| , and is illustrated in the inset 
of Fig. O We note that the results presented here do not 
depend significantly on the parameters which describe 
the electrostatic walls. The weak disorder potential used 
in the simulations is modeled as a Gaussian random field 
with zero mean and standard deviation of the amplitude 
of the disorder potential Vq = 2%Ei?, and a Gaussian 
correlation function (V(r)V(rO) = Fo^e-l^-^'l'/^'- , with 
correlation length I cor = 180 nm. As these parameters 
are not readily accessible in the experiment, they were 
chosen from a realistic range of values^. The mean free 
path in this model potential is even larger than the mea- 
sured mean free path, which is limited by other scatter- 
ing processes, e.g. scattering by charged crystal defects. 
These are, however only relevant on length scales larger 
than our system size. Although our model system is even 
"deeper" in the ballistic regime than suggested by the 
experimentally evaluated mean free path, the weak dis- 
order nevertheless has a pronounced effect on the peak 
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FIG. 6. Quantification of the influence of disorder, (a) The 
average distance an electron travels until a caustic forms as 
a function of the disorder parameters for B — 0^ B — Bq 
and B — 2Bo. The dotted lines mark the parameters of the 
random potential chosen for the simulations in Fig. \5\ (b) The 
numerically calculated inverse participation ratio (IPR) of the 
flow density along the bottom boundary between x = 0.5a and 
X = 1.5a, where a is the distance from QPC to the corner, as 
a function of the standard deviation of the amplitude of the 
disorder potential under coUimation and focusing conditions 
(B = Bo and B = 2Bo), averaged over 200 realizations of the 
random potential. The correlation length is identical to the 
one for the calculations in Fig. [S] and the IPR is normalized 
by the IPR of the clean system. The rising IPR oX B = Bq 
indicates that increasing disorder produces more peaks in the 
flow density, while the focusing peak B — 2Bo is broadened 
by the disorder. 



structure. Figure [5^ shows calculated resistance traces 
for different realizations of disorder potential with the 
same statistical parameters. We point out that the re- 
sults are not sensitive to small variations in the param- 
eters chosen for the random potential. Figures [Sb and 
[Sb compare the experimental traces from Figs.|4]3 andHJ) 
with the calculated resistances. Finding potential disor- 
der landscapes which produce such excellent agreement 
with experiment is to some extent accidental. The main 
purpose of the simulations is to demonstrate that differ- 
ent disorder landscapes characterized by the same statis- 
tical parameters can indeed either lead to a pronounced 
splitting of collimation related peaks as one would ex- 
pect from the behavior of the flow density at the sample 
boundary in the bottom left panel of Fig.[3li, or not. The 
focusing features at even values of B / Bq do not show 
any splittings for the chosen parameters of the disorder. 
With decreasing correlation length or increasing Vo how- 
ever, we anticipate that also magnetic focusing features 
are affected by disorder induced branching. To assess the 
influence of the disorder as a function of the correlation 
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length and Vq, it is instructive to calculate the average 
distance an electron travels until a caustic forms, IcausUc 
We have obtained an analytical expression for this quan- 
tity, which is given by 

Icaustic = TT^K (Ai [-k2b2]2 + [-^^B^]') 

where Ai and Bi are Airy functions of the first and sec- 
ond kind and where is a function of the disorder po- 
tential (see Appendix A). Our result is plotted in Fig. [6^ 
for 5 = 0, 5o and 2Bq. For the latter field, magnetic 
focusing causes a caustic in our corner geometry after 
electrons have traveled on average a distance 7ra/2. In- 
deed, Icaustic saturates to this value in the limit of large 
correlation lengths and small Vq. For both values of the 
magnetic field we observe that for the parameters chosen 
in the simulations (indicated by dashed vertical lines) the 
mean distance to the first caustic starts to deviate from 
the case of magnetic focusing without disorder. For zero 

magnetic field Icaustic scales like /^or/V^^^^'^*^*^ Hence, 
branching can occur on much shorter length scales than 
the mean free path, which scales as I^q^/V^^ Our an- 
alytical calculation shows that branching influences the 
transmission properties of our device. From Fig. [6^ it 
is however not apparent that disorder has a different 
impact on collimation [B = Bq) than on focusing fea- 
tures {B = 2^0 )• To assess the impact of branching 
more quantitatively, we study the peakedness of the flow 
density along the bottom boundary of the corner de- 
vice. Examples of the flow density I{x) are displayed 
in Fig. [3b and d. As evident from this figure, the dis- 
order potential induces multiple peaks near the bottom 
QPC dX B = Bq^ while the peakedness of the flow den- 
sity for B = 2Bq resembles that of the flow density in 
the absence of disorder. To capture the peakedness of 
the curves in a single quantity we use the inverse par- 
ticipation ratio IPR[I]=/^''^' dxP{x)/{j^l dxl{x)f. This 
inverse participation ratio in the vicinity of the bottom 
QPC is plotted as a function of the standard deviation 
of the disorder amplitude in Fig. [6)3. For the collima- 
tion condition B = Bq, IPR[I] rapidly rises to higher 
values. It reflects the appearance of additional peaks. 
Under focusing conditions, IPR[I] drops indicating that 
the original peak mainly broadens. This confirms our 
experimental observations. 

V. CONCLUSION AND ACKNOWLEDGMENT 

In conclusion, the weak disorder potential, which is 
present in any real two dimensional electron gas, causes 
a pronounced modification of collimation features due to 
branching. Our findings indicate that branching needs to 
be taken into account when interpreting transport data of 
mesoscopic devices even for state-of-the art heterostruc- 
tures, for which the mean free path is an order of mag- 
nitude larger than the device size. 
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Appendix A: Mean distance to the first caustic in a 
random potential and magnetic field 

In order to determine the location of caustics, we con- 
sider an equation for the curvature of the action function 
/S, which is obtained from the Hamilton- Jacobi-equation 
(HJE). Points along a trajectory where the curvature 
diverges indicate the position of a caustic. In a con- 
stant magnetic field 5, the electron trajectories are circu- 
lar. Considering small deviations from the circular paths 
in polar coordinates {r, (j)} allows a quasi-2D treatment, 
similar to the one of Fig. 2, in which time is identified 
with the angular variable as follows. The HJE in the 
symmetric gauge with vector potential A = ^r^ is given 

by 

dtS + ^{drSf + ^r^B^ + V{r{t)) = 0, 

where we have identified t with j . Taking two derivatives 
with respect to r , and evaluating the equation for the 
curvature u = drrSir) along the trajectories, we obtain 
the following equation for u: 

j^u^u^ ^B'^ ^ drrV{r) = 0. 

For weak random potentials, we can approximate the 
correlation function of the random potential as c(r — 

r',(t)-(t)') = (y(r,0)U(r>O) =_i(£.^0O^(^ - ^0- 
Thus, extending results from Refs. [27.30.31], we derive a 
Fokker-Planck-equation for the probability density p(t, u) 

where D = -^c{x,y)\y=odx. For a Gaussian corre- 
lation function 

we obtain D = 12^/7TVQl~^^. To obtain an equation for 
the onset of the branching, we now derive an expres- 
sion for the mean time it takes for a caustic to develop 
along a trajectory. This can be done treating the prob- 
lem as a mean first passage time problem [26||32|. From 
this expression, one can then easily derive the mean dis- 
tance traveled along a trajectory until a caustic is hit. 
Of greatest importance for the experiment is the point 
source with initial condition uq = oo. The corresponding 
mean distance to the first caustic Icaustic is then calcu- 
lated [26*] to be 

lcaust^c = ^ (aI [-^^B^] ' + Bi [-K^B^] 

where Ai and Bi are Airy functions of the first and second 
kind m, and where = 
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